Scaling Properties of Growing Noninfinitesimal Perturbations in Space-Time Chaos 
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We study the spatio-temporal dynamics of random spatially distributed noninfinitesimal pertur- 
bations in one-dimensional chaotic extended systems. We find that an initial perturbation of finite 
size eo grows in time obeying the tangent space dynamic equations (Lyapunov vectors) up to a 
characteristic time ix(eo) ~ b — (1/Xmax) ln(eo), where \max is the largest Lyapunov exponent 
and b is a constant. For times t < t x perturbations exhibit spatial correlations up to a typical 
distance £ ~ t z . For times larger than t x finite perturbations are no longer described by tangent 
space equations, memory of spatial correlations is progressively destroyed and perturbations become 
spatiotemporal white noise. We are able to explain these results by mapping the problem to the 
Kardar-Parisi-Zhang universality class of surface growth. 
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I. INTRODUCTION 

A standard tool for studying chaotic behavior in dy- 
namical systems is the computation of the characteristic 
Lyapunov exponents, which measure the typical expo- 
nential growth rate of an infinitesimal disturbance 0,13 • 
The characteristic Lyapunov exponents in extended sys- 
tems are defined in a similar way as their low-dimensional 
counterpart and can be calculated from the linearization 
of the equations of motion 0, ■ The main point is that 
the growth of an infinitesimal perturbation is described 
by the linear equations for the tangent space, the so- 
called Lyapunov vectors (see below). However, for many 
practical purposes the dynamics of infinitesimal pertur- 
bations may be irrelevant as indicators of, for instance, 
the predictability time. Indeed, in many realistic situ- 
ations the error in the initial conditions is finite. The 
important fact is that the evolution of finite errors is not 
confined to the tangent space, as defined by the growth 
of linearized perturbations, but is controlled by the com- 
plete non-linear dynamics. A good example with impor- 
tant practical application occurs in weather forecasting, 
where one deals with the whole Earth's atmosphere- an 
extremely high dimensional system in which initial condi- 
tions can be determined only with limited accuracy. The 
effects of finite perturbations have recently been studied 
in the context of fully developed turbulence In or- 
der to deal with realistic perturbations, the concept of 
finite size Lyapunov exponents has been found to be use- 
ful to analyze predictability in high-dimensional chaotic 
systems 0,0- 

In this paper we study the dynamics of random spa- 
tially distributed finite-size errors in chaotic extended 
systems and focus on their propagation dynamics. We 



argue that, after a suitable transformation of variables, 
the dynamics of finite perturbations can be interpreted as 
a kinetic roughening process in the Kardar-Parisi-Zhang 
(KPZ) universality class [l0|. We find that, due to the 
finiteness of the initial error, there is a characteristic time 
scale t x (e ) ~ b - (l/X ma x) M e o), where X rnax is the 
largest Lyapunov exponent, eo is a measure of the initial 
size of the perturbation, and b is a constant. For times 
t < t x the dynamic evolution of a finite perturbation is 
governed by the Lyapunov vector. In this regime, finite 
perturbations become spatially correlated up to a typical 
length scale £ ~ t z x , where z is the dynamic exponent of 
the KPZ problem (z — 3/2 for one dimensional systems). 
However, for times t > t x any finite perturbation leaves 
the tangent space and is no longer described by the Lya- 
punov vectors. In this late regime, memory of spatial 
correlations is progressively destroyed and perturbations 
actually become white noise in space and time. Our ap- 
proach provides new tools for studying chaotic extended 
systems by allowing to fully describe the spread of corre- 
lations, to estimate the spatial extend of correlations of 
the chaotic field and to measure the effective number of 
degrees of freedom. 

We exemplify our results by means of numerical simu- 
lations of coupled map lattices in one dimension, which 
are simple model systems exhibiting space-time chaos 
and convenient as far as the computing time is concerned. 
We consider a coupled map array consisting of L chaotic 
oscillators given by 



u(x,t+l) = uf(u(x + l,t)) + 
vf(u(x-l,t)) + (l-2v)f(u(x,i)) 



(1) 
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where x — 1, 2, . . . , L, f(u) is a chaotic map, v is the 
coupling constant, and periodic boundary conditions are 
imposed. We have fixed the coupling to v = 1/3 in all 
the simulations presented in this paper. We have carried 
out simulations for two different choices of the map, the 
chaotic logistic map f(u) = 4u(l — u), < u < 1 and the 
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tent map f(u) = 1 — 2\u — 1/2|, < u < 1. For the sake 
of brevity, all the results we present below correspond to 
coupled logistic maps, but similar results were obtained 
for the tent map. 

The dynamics of infinitesimal perturbations of a tur- 
bulent state u°(x, t) can be studied by linearizing the evo- 
lution equation Eq. Q). This leads to the tangent space 
equations or Lyapunov vector 5u(x, t) evolution equation 

Su(x, t + l)=v f[u{x + 1, t)] Su(x + l,t) + 
vf[u{x- l,t)]5u(x-l,t) + 
(1 - 2 v) f'[u(x, t)} Su(x, t) + 0[(5u) 2 } (2) 

where f'[u(x,t)] — df (y) / dy\ u ( x , t ) ■ This implies to solve 
simultaneously the field u(x,t) evolution equation 
The Lyapunov vector is defined in the linear approxi- 
mation, when higher order corrections 0[(<5it) 2 ] are ne- 
glected. 

The analysis of infinitesimal perturbations allows to 
compute several indicators characterizing the chaotic sys- 
tem, including the whole spectrum of Lyapunov expo- 
nents and investigate how this depends on system size 
0L1E1 However, as explained above, there are indeed 
many situations where the Lyapunov analysis has no rel- 
evance due, for instance, to the finite nature of the initial 
errors. 



II. SCALING OF FINITE-SIZE 
PERTURBATIONS 

Let us now consider the evolution of random finite-size 
perturbed trajectories in our model system Q). Given 
an initial condition m°(x,0), the solution vP(x,t) is uni- 
vocally determined by computing Eq. Q for a number t 
of time steps. This is our reference trajectory and we are 
interested in the evolution of finite perturbations around 
that reference solution. For reasons that will become 
clear below, we find convenient to introduce now what 
we call the amplitude factor e(t) as the spatial geometri- 
cal mean value of a perturbation: 

L 

e(t) = J[\5u(x,t)\ 1 ' L . (3) 

x=l 

As we shall see below, the amplitude factor turns out to 
be a very important quantity which contains the infor- 
mation about the dominant exponential growth rate. 

Since we are interested here in the propagation of real 
(non-infinitesimal) errors, we should avoid linearization 
of Eq. Q. Instead, we compute the trajectories gen- 
erated by iterating Q for an ensemble of initial con- 
ditions u(x,0) = u°(x,0) + Su(x,0), where Su(x,0) is 
a random finite perturbation of initial amplitude e(t — 
0) = e . For each iteration of the lattice the dif- 
ference Su(x, t) — u(x, t) — t) between the reference 
trajectory and every one of the perturbed solutions is cal- 
culated. Although the disturbances are initially random 
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FIG. 1: Numerical results for the propagation of finite- 
size errors in coupled logistic maps. Upper panel shows 
the ensemble averaged amplitude factor (e(t)) vs. time 
for perturbations starting with initial amplitudes of eo = 
10" 3 , 10" 5 , 10" 7 , 10~ 9 , 10~ n and 10" 13 (from top to bottom) 
in ID lattices of L — 1024 sites. Results were averaged over 
600 different initial conditions. Lower panel shows the vari- 
ance a 2 (t) for the same initial perturbations as before and 
eo decreasing from 10 -3 (leftmost curve) to 10 -13 (rightmost 
curve) . 



and uncorrelated in space, as time goes by, they grow 
and get spatially correlated. The statistical fluctuations 
of errors can be characterized by studying the ensemble 
of finite perturbations {5u n (x,t)}^ =1 , which correspond 
to N independent realizations of the initial perturbation. 

In Fig. 1 (upper panel) we plot ln(e(i)) vs. time for 
different values of the initial perturbation amplitude eo, 
where (• • ■ ) stands for average over realizations of the ini- 
tial perturbation. One can immediately see that there ex- 
ists a characteristic time scale t x (eo) such that for times 
t < i x (eo) the average amplitude factor grows exponen- 
tially in time e(t) « eo exp(Ai), where A = 0.343 ± 0.005. 
We demonstrate below that A indeed corresponds to the 
maximal Lyapunov exponent. For longer times, t > 
tx (eo) the amplitude factor saturates to a constant value. 
Both, the saturation constant and the maximal Lyapunov 
exponent, are independent of the initial perturbation size 
eo- However, the saturation times i x (eo) increase as the 
size of the initial perturbation eo becomes smaller. The 
characteristic time scale i x corresponds to the crossover 
time at which the dynamics of a finite size perturbation 
depart from the linear evolution (i.e. the Lyapunov vec- 
tors) given by Eq. This crossover occurs because 
Lyapunov vectors describe only the behavior of strictly 
infinitesimal perturbations. One can estimate t x roughly 
as the time at which e(t) reaches some 5, so that the 
higher order terms 0[(<5u) 2 ] cannot be neglected in the 
evolution equation @. This crossover takes place at a 
typical time t x (e ) - (1/A)ln(<5) - (1/A)ln(e ). There- 
fore, for times t > i x (eo) nonlinear corrections, due to 
finiteness of the initial perturbation, come into play and 
drive errors out of the tangent space. From then on, 
the linear approximation cannot describe the evolution 
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of errors. One then expects that t x (eo) oo as e ^ 0. 

Besides exponential growth, spatial correlations are 
dynamically generated during the evolution of the per- 
turbation. Correlations contain information about the 
sub-leading Lyapunov exponents and thus also contribute 
to the perturbation size growth. The important role of 
correlations can be better realized after subtraction of 
the dominant exponential growth component given by 
e(t). We find that a very useful indicator is given by the 
reduced perturbations 8r(x, t) that we define as 



where the dominant exponential growth is globally re- 
moved and one is left with the effect of correlations. For 
random initial errors, i.e. random spatially distributed 
initial conditions, the average (Sr(x, t)) vanishes, where 
(■••} stands for average over realizations of the initial 
perturbation and the over bar is a spatial average. Thus, 
statistical fluctuations of the reduced perturbations are 
measured by the variance <r 2 (t) = (Sr(x, t) 2 ). As we shall 
see below, a(t) gives information about the growth of 
perturbations due solely to correlations. 

In Fig. 1 (lower panel) we show our numerical results 
for the variance of the reduced perturbations, Eq.Q). 
The variance <r 2 (t) grows as <J 2 (t) ~ exp(i 2/3 ) for times 
t < t x , i.e. before saturation occurs due to finiteness of 
the initial disturbance. This scaling behavior is better 
seen in Fig. 2 where we plot in log-log scale ln(er 2 ) vs. 
time and the exponent f3 — 0.30 ± 0.05 (see Fig. 2). As 
before, this rapid growth occurs for times t < t y when 
the dynamics of disturbances are well described by the 
Lyapunov vectors. One can determine the crossover time 
t x from Fig. 1, either by measuring the saturation times 
of e(t) (upper panel) or by measuring the time shift of 
the maxima of a 2 (t) (lower panel) . Results are shown in 
the inset of Fig. 2, where a straight line fit of the data 
leads to £ x (eo) = b— (l/A)ln(eo) with b a constant. The 
slope of the fit 1/A = 2.86 ±0.05 is in excellent agreement 
with our previous determination of A w 0.34 from the ex- 
ponential growth of the amplitude factor e(t) ~ exp(Ai). 

The magnitude and extent of spatial correlations can 
be measur ed by using the site -si te correla tion function 
G(l,t) = (Sr(x ,t)Sr{x + I, t))/{Sr(x ,t) 2 ). The mag- 
nitude of spatial correlations increases in time G{l,ti) < 
G(l, £2) if t\ < ti for times tj., ta < t x . However, as shown 
in Fig. 3, correlations become progressively smaller for 
times t > t x . This again demonstrates that the effect of 
having initially finite perturbations is to introduce a char- 
acteristic time f x (eo) marking the typical time it takes 
for the system to depart from tangent space (with the 
building-up of correlations) to non-linear evolution (un- 
correlated errors). For times larger than t x the pertur- 
bation spatial correlation are progressively destroyed, as 
can be seen from the decay of the correlations G(l,t) 
in Fig. 3. Further Fourier spectrum analysis shows 
that errors actually become undistinguishable from white 
noise. 




FIG. 2: The same numerical results for the variance <r 2 (t) as 
those in Fig. 1 are plotted in the main panel in log-log scale 
to show the power-law behavior lna 2 (t) ~ t 213 . The straight 
line has a slope 0.60 and is plotted to guide the eye. The 
inset shows the mean-square fit of the characteristic time scale 
ix(eo) obtained from the saturation times of the amplitude 
factor e(t). The slope of the straight line is 1/A = 2.86(5). 



III. KINETIC ROUGHENING PICTURE 

In order to make a theoretical interpretation of our 
numerical results, in particular the existence of scaling 
behavior, we can resort to the Hopf-Cole transformation, 
first proposed by Pikovsky and Politi for the actual Lya- 
punov vectors jl2l ll3T | , that allows to link spreading of 
errors to non-equilibrium surface growth. We find that 
the dynamics of finite perturbations can also be seen as 
a kinetic roughening processes of the surface defined by 
h(x,t) = hi\5u(x 1 t)\. This is a very useful transforma- 
tion that allows us to borrow some well-known results 
in the field of nonequilibrium surface growth. Indeed, 
Politi and Pikovsky [T2I Il3l ] have shown that errors in 
many extended systems lead to surface growth process 
in the universality class of KPZ [Io|. 

Let us first consider the amplitude factor defined in 
Eq. @ and show that it is related to the average sur- 
face velocity. We can write |<5u(a;, i)| = exp[h(x,i)] and 

thus e(t) = exp[(l/L)Y%Zi h(x,t)] = exp[7i(*)]. Wc 
then obtain that the amplitude factor must grow as 
e(t) ~ eo exp(A maa; t), where X m ax is the largest Lyapunov 
exponent and corresponds in this mapping to the surface 
velocity, h(t) = X max t + h(Q) 0,0- Direct measures of 
the average surface velocity (not shown) are in excellent 
agreement with our above discussed estimations of A. As 
an independent check we have also measured the largest 
Lyapunov exponent by standard techniques [hH Il5j . 

Also the time behavior of the variance of the re- 
duced perturbations shown in Fig. 2 can be linked to 
surface scaling properties. From the surface definition 
h(x,t) = In \5u(x, t)\ we have a 2 (t,L) = (exp[2y(x, £)]), 
where y(x,t) = h(x,t) — h(t). We can now proceed and 
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FIG. 3: Time behavior of the correlation function for lattice 
sites separated at different distances I = 8 (o), I = 16 (□), 
I = 64 (<), and I = 185 (>) in an array of L — 1024 coupled 
maps. 



Lyapunov exponents exists and whose number increases 
with the system size. 

The ansatz of a KPZ behavior for the surface h then 
implies that the width should scale as W(t, L) ~ t@ for 
times t < t s (L) and saturates to a size dependent value, 
W(t, L) ~ L a for t > t s (L), where f3 and a are the growth 
and roughness exponent, respectively, that are known to 
take the values /3 = 1/3 and a = 1/2 for the KPZ univer- 
sality class in one dimension 0, ^} . This is consistent 
with the approximate exp(i 2 / 3 ) scaling observed in Fig. 
2, baring in mind the limitation of the approximation 
©• It is worth mentioning that all the numerical re- 
sults presented here are obtained for system sizes such 
that t x t s , so that the existence of t x could be clearly 
observed. 



calculate a 2 explicitly by considering the cumulant gen- 
erating function of the stochastic variable y given by 
$(s) = ln(exp(zsy)} |Tfjj . where i is the imaginary unit. 
This generating function can be expanded in a power se- 
ries of the argument s and then evaluated at s = — 2i to 
obtain an exact expression for the variance a 2 : 

^,L)=exp[f>^], (5) 

r=l 

where {{y r )) are the cumulants of y [Tf^ . In general, 
no closed simple formula exists for the cumulant of or- 
der r of a given variable y, but the first cumulants 
are_ given by ({y))_= (y), ((y 2 )} = (y 2 ) - (y) 2 , and 
((y 3 )) = (V 3 ) - 2(y 2 )(y) + 2(y) 3 , etc 0. In our case 
(y) = and the second cumulant ((y 2 )) — W 2 (t,L) cor- 
responds to the surface width W 2 (t, L) — ([h(x, t) — h] 2 ). 
To leading order we have 

a 2 {t,L) w exv[2W 2 (t,L)]. (6) 

Higher order cumulants are also finite, but their contri- 
bution to the series gets smaller as r increases. We have 
checked numerically that including the 3rd and 4rd cumu- 
lants shifts the exponent in less that a 10%. Therefore, 
within this somewhat crude approximation, Eq.JfjJ), we 
expect to have an estimation of the effective scaling ex- 
ponent of a 2 ~ exp(t 2 ^) with less than a 10% error bar. 
A clarification is now in order. For systems in which 
the largest Lyapunov exponent A — > one would expect 
to have larger crossover times tx, so that higher order 
cumulants may eventually become relevant, yielding cor- 
rections to Eq. ©. In this sense the scaling picture has 
to be considered as approximate. However, the scaling 
approach is consistent for hyper-chaotic systems, as those 
we are interested here, in which not only the leading Lya- 
punov exponent is finite, but a series of finite positive 



IV. CONCLUSIONS 

We have studied the spatiotemporal dynamics of finite- 
size perturbations in extended systems exhibiting chaos. 
We have introduced what we call the amplitude factor 
e(i) in Eq.® and showed that this quantity is directly 
related to the maximal Lyapunov exponent and exhibits 
nice scaling properties. The amplitude factor allows us to 
remove in a simple and self-consistent way the dominant 
growth component of the perturbation at every time step, 
so that spatial correlations can be determined more eas- 
ily. At variance with infinitesimal perturbations, finite- 
size errors lead to a characteristic time scale t x {co) sig- 
naling the departure of the finite perturbation from the 
linear approach. For times t > t x the evolution is nonlin- 
ear and spatial correlations are destroyed. We have found 
that the perturbation variance grows as a(t) ~ exp(£ 2//3 ), 
once the contribution of the largest growth rate is re- 
moved. We have explained this behavior by representing 
the error growth as a kinetically rough ID KPZ surface. 

The interpretation of the exponentially growing errors 
Su(x, t) as a roughening surface h(x, t) = In \5u(x, t)\ has 
many advantages since it allows to translate the problem 
of analyzing deterministic chaotic fluctuations into the 
simpler framework of kinetic roughening of the associ- 
ated surface h(x,t). The scale-invariant character of the 
surface fluctuations leads to simple power-law behavior 
of the relevant quantities that can be characterized by a 
few critical exponents. As a byproduct this leads in a 
rather simple way to an estimation of the spatial extent 
of correlations at time t as £ ~ t z , where z = 3/2 is the 
dynamic exponent of the KPZ problem. Further results 
of practical application, including the relevance for the 
predictability problem of the propagation of finite-size 
errors and bred vectors, as used in the context of weather 
forecasting will be published elsewhere [Tl| . 
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